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The width of three-body single-pion decay process d* + NNrx°* is calculated by using the d* 
wave function obtained from our chiral SU(3) constituent quark model calculation. The effect of the 
dynamical structure on the width of d* is taken into account in both the single AA channel and cou- 
pled AA+CC two-channel approximations. Our numerical result shows that in the coupled-channel 
approximation, namely, the hidden-color configuration being considered, the obtained partial decay 
width of d* + NWNz is about several hundred KeV, while in the single AA channel it is just about 
2~ 3 MeV. We, therefore, conclude that the partial width in the single-pion decay process of d* 
is much smaller than the widths in its double-pion decay processes. Our prediction may provide a 
criterion for judging different interpretations of the d* structure, as different pictures for the d* may 
result quite different partial decay width. 
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In recent years, CELSIUS/WASA and WASA@COSY Collaborations (a) [py have reported a clear evidence of a 
resonance-like structure in double pionic fusion channels pn > dr°n° and pn + drr when they studied the ABC 
effect [3] and when they treated the neutron-proton scattering data with newly measured analyzing power A,. Since 
the observed structure cannot be simply understood by either the intermediate Roper excitation contribution or the 
t-channel AA process, they proposed an assumption of existing a d* resonance whose quantum number, mass, and 
width are (JP) = 0(3+), M ~ 2370 MeV and T = 70 MeV [il [Qj (see also their recent paper [4], the averaged mass and 
width are M ~ 2375 MeV and I = 75 MeV, respectively). Due to its baryon number being 2, it would be treated as a 
dibaryon, and could be explained by either ” an exotic compact particle or a hadronic molecule” [B]. Moreover, accord- 
ing to the experimental data, the mass of the d* is about 80 MeV below the AA threshold and about 70 MeV above 
the ATN threshold, therefore, the threshold (or cusp) effect is expected not to be so significant as that in the XYZ 
study (see the review of XYZ particles (6]).Thus, understanding the internal structure of d* would be of great interest. 


Actually, the existence of such a non-trivial six-quark configuration with (JP) = 0(3+) (called d* lately) has 
triggered a great attention and has intensively been studied in the literature since Dyson’s pioneer work . It 
should specially be mentioned that one of those calculations reported in 1999 studied the binding behavior of the 
37 dibaryon system by taking into account a AA channel and a hidden-color channel (denoted by CC hereafter) 
simultaneously m. In that paper, the binding energy was predicted to be about 40 — 80 MeV which is quite close 
to the recent observation, and the importance of the CC channel was particularly emphasized. Unfortunately, the 
width of the state was not calculated. 


After the discovery of d*, there are mainly three types of explanation for its nature. Based on the SU(2) quark 
model, Ref. proposed a AA resonance structure and performed a multi-channel scattering calculation. They 
obtained a binding energy of about 71 MeV (namely Ma» = 2393 MeV) and a width of about 150 MeV which is 
apparently much larger than the observation. On the other hand, Ref. Ad studied a three-body system of AN7 and 
found a resonance pole with a mass of 2363+20 MeV and a width of 65+17 MeV. An important view point, claimed by 
Bashkanov, Brodsky and Clement in 2013, is that a dominant hidden-color structure (or six-quark configuration) 
of d* is necessary for understanding the compact structure of d*. Sooner after, following our previous prediction í 
Huang and his collaborators made an explicit dynamical calculation by using a chiral SU(3) constituent quark 
model Edis in the framework of the Resonating Group Method (RGM), and showed that the d* state has a 
mass of 2380 — 2414 MeV, which agrees with COSY’s observation, and does have an explicit “CC” configuration 
of about 66 — 68% in its wave function ag. Based on the obtained wave functions of d* and deuteron, Dong and 
his collaborators calculated the partial decay widths of the Golden” decay channel d* > d + 27° (ntr) ba) as 
well as the widths of its four-body decay d* — pnx°7r° and d* + pnrtr [21]. The results of the two papers 
showed that inclusion of the CC configuration inside d* would make the calculated widths suppressed greatly and 
the resultant partial widths for all the d* —> dr°r°, d* 3 drr”, d* > pnr°r®, and d* > pnatn~ decay processes 
are well consistent with the experimental data. As a consequence, the total width of d* is about 72 MeV, which 
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is fairly close to the observation. All these outcomes support that d* is probably a six-quark dominated exotic 
state due to its large CC component. The general review on the dibaryon studies can be found in Ref. by Clement. 


Recently, the questions about how large the decay width of the single pion decay mode of d* is and whether such 
a decay process can be observed have been discussed. Up to now, the d* — N Nr decay process has not been found 
in the data analysis 23h, but_a sizable cross section of this process was predicted by using a ANa model where the 
rms of d* is about 1.5 fm Rd. This contradictory information encourages us to calculate the partial decay width 
of this process in the same way employed in calculating partial decay widths of the dma and NNrr processes and 
with the same d* wave function obtained in our AA + CC model, with which the resultant rms of d* is only about 
0.8 fm aig, and all the calculated partial decay widths of the drr and NNrr processes and the total width of d* 
are consistent with the data quite well [20, Bi. Therefore, the obtained partial width might be used to distinguish 
the models for the structure of d*. 


Similar to our previous work 2d, Bi, we employ the phenomenological effective Hamiltonian for the pseudo-scalar 
interaction among quark, pion, and quark in the non-relativistic approximation 


= — Iuar F. je . 
Haqr meg (2T)3/2 V2wr o kinT Q, (1) 


where gqqr is the coupling constant, @ stands for the m meson field, wy and kr are the energy and three-momentum 
of the m meson, respectively, and o(r) represents the spin (isospin) operator of a single quark. The wave func- 
tions of the nucleon and A(1232) resonance in the conventional constituent quark model can be found in 2o). 
The experimental data for the decay width of A —> aN is ~ 117MeV, and the theoretical calculation gives 
Pasan = sk? (Jagr lo) En /Ma, where Ey = y Mẹ, + py is the energy of the outgoing nucleon, ky ~ 0.229 GeV is 
the three-momenta of pion, and J, denotes the spatial overlap integral of the internal wave functions of the nucleon 
and the A resonance, we can extract the coupling constant gqgqr (the details can be found in Ref. [20}). Defining 
G = gqqr Io, the obtained G value is about 5.41 GeV~!. 


As mentioned in Refs. [19, Bo, our model wave function is obtained by dynamically solving the bound-state RGM 
equation of the six quark system in the framework of the extended chiral SU (3) quark model, where the one-gluon- 
exchange and Goldstone Boson exchange interactions between quarks are explicitly considered. Then, by projecting 
the obtained wave function onto the inner cluster wave functions of the AA and CC channels, the wave function of 
d* can be abbreviated to a form of 


Was = [ ba(&, i) palé E) xaa(R) Can + b0(61,€) pelë, 6) xee(R)) Coe |(sn=(30); (2) 


where ġa, and ¢c denote the inner cluster wave functions of A and C (color-octet particle) in the coordinate space, 
xaa and xycc represent the channel wave functions in the AA and CC channels (in the single AA channel case, 
the CC component is absent), and aa and cc stand for the spin-isospin wave functions in the hadronic degrees 
of freedom in the AA and CC channels, respectively [1 9}. It should be specially mentioned that in such a d* wave 
function, two channel wave functions are orthogonal to each other and contain all the totally anti-symmetrization 
effects implicitly Eg). 


In terms of the obtained d* wave function eq. (2) (its explicit forms have been plotted in Ref. [20}) and, we are able 
to calculate the width of the three-body single-pion decay d* = NN. The partial width reads 


2 


1 oe 
Pann = 5 | Prd’ 2n)s(AE)| MOP > (3) 


2 
| stands for the squared transition matrix element with a sum over the polarizations of 


where [MG Ps) 
the final three body states and an average of the ones of the initial state d*, the factor of 2! is due to the 
property of the identical particle in the final NN system, and (AE) denotes the energy conservation with 
AE = Ma — wrk) — En(p1) — En (—pi — k) where w,(k) and Eyn represent the energies of the pion and nucleon, 
respectively. 


The transition matrix contains contributions from 12 Feynman diagrams, where the AA component in the 
d* wave function is responsible for the decay, the pion-exchange is considered for those sub-leading effects, and 
the intermediate nucleon state is taken into account only. Among these diagrams, 6 of them where the outgo- 
ing pion is emitted from A» are drawn in Fig.1, and they are depicted according to the time-order perturbation theory. 
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FIG. 1: Six possible ways to emit pion only from the AA component of d* in the d* — N Nz decay process. The outgoing 
pion with momenta k is emitted from A2. The other six sub-diagrams with pion emitted from Aj are similar, and then are 
not shown here for reducing the size of the figure. 


Computing transition matrix elements for all the diagrams in Fig. 1 is straight forward (refer to the discussions in 
Ref. ). For example, the matrix element for Fig. 1(a) can be written as 


Vax (q ) 
Mo = fe —— eit __53 r+ 1+ k— = 
ENNE lwr, Vor T)’ (pws RNI På Pas) 


xMarlke) No) =N (p2) Dat Ma, —+m(kn)N(p1)DaiMas+m(k)N (ph) (4) 


where Wy» represents the d* wave function in the momentum space which can be obtained by Fourier transforming 
the d* wave functions in the coordinate space in both the single AA channel and coupled AA + CC channel approx- 
imations (see details in Refs. ) Ma(ke)N (ph) >N (p2) May sa(ke)N(p1)> and Ma,-+x(k)N (pi) denote the transitions 
of t(kp)N(p5) > N(p2), Ar > 7(kge)N (pı), and A2 > x(k) N(p5), respectively, and the non-relativistic energy 
propagators are 


Dp = ——— SS 
Mar — w(k) — w(ke) — En (pi) — En (py ) 


e E (5) 
Ma» — w(k) — Ea, (7) — En (Dz ) 


Then the explicit form of the matrix element (Fig. 1(a)) for the case where the spin of the final two-nucleon is zero 
(s12 = 0) can be expressed as 


GannG2 at 2V2 s [poses qWax(q)kz, xy 3 
k) = ——— E — kY o(k Yom, ( k SON 
Ma- o(P1; ) (27r)8/2u%, 27 1,0( ) eae D 2 i E) 5 C10, 2m (6) 


In the above equation, ke = q — pı stands for the three-momenta of the exchanged pion, Gann & 4V2G and 


Gana ~ “yes (with gvnz ~ 13.6), and kp = @— pı. Moreover, for the case of s12 = 1, we have 


_ GrunGeant 2V2 
(2r) V2wk 27 


1 jim 2 yl 
s12 3Max Meig 3Max 
x [y ia Cimen 52m) g5 3O10,.2me1, Comen ™ è (7) 


kY, o(k) Pava la)ki ye Yom, (ke) 


Mo, k = 
(prk) 2wke Dai Daf 
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In general, the final state interaction (FSI) between two outgoing nucleons should be considered. However, being 
aware of the fact that the quantum numbers of d* are I(J”) = 0(3*), the maximal spin of two nucleons is 1, and 
parity P and total angular momentum J of the decaying system should be conserved, either the orbital angular 
momentum between the outgoing pion and nucleon is at least equal to 3, or the orbital angular momentum between 
two outgoing nucleons at least equals to 2. In the former case, the decay cross section would be greatly suppressed 
by the higher partial wave. And in the latter case, the FSI effect could be ignored because in the low energy region, 
the dominant contribution comes from the S-wave and P-wave, and the contribution from the higher partial wave 
can be ignored (26). Therefore, in this calculation, we assume the enhancement factor from FSI is close to 1, and 
consequently would not be considered. 


Table 1, The calculated decay width of the d* — pnr? process and the widths contributed individually 
from the (a)-, (b)-, (c)-, (d)-, (e)-, and (f)-type diagrams (in units of MeV). 


Gwe ido) | oO) oO] @ | oO) ] am ora) 


[Case 
One ch. (AA only) | 2276 || 0.550 | 0306 0.0963 
Two chs. (AA+CC) || 0.670 || 0.154 | 0.0884 | 0.0789 | 0.0279 | 0.0687 | 0.0847 0.503 


The numerical results for the decay width in the d* — pnn? process and the widths contributed individually from 
the (a)-, (b)-, (c)-, (d)-, (e)-, and (f)-type diagrams are tabulated in Table 1, respectively. Contributions from the 
diagrams where the outgoing pion being emitted from Aj are also taken into account. The total decay width of 
the d* — NWNr7 process is the sum of the contributions from all types of diagrams coherently. From table 1, one 
sees that the ratio of the decay width with coherent sum to that with incoherent sum is about 1.37 in the one 
channel (AA channel only) case and about 1.33 in the two channel (coupled AA+CC channels) case, respectively, 
which shows a sizeable coherent effect. The most important issue from this table is that the decay width of d* to 
pnr” is smaller than 3 MeV in the one channel case and about 670 KeV in the two channel case. Apparently, the 
width in the two channel case is much smaller than that in the one channel case. Since in the two channel case, the 
contribution comes from the AA component of d* only, which is about 31.5% of the whole d* wave function, but 
in the one channel calculation, the contribution comes from the whole AA wave function, therefore, this outcome 
is understandable. The second observation is that although the framework and method in this calculation are the 
same as that used in the d* width calculations in the double-pion decay processes before, the obtained decay widths 
for the single-pion decay mode are remarkably smaller than those given in our previous calculations (20, and the 
experimental data for the double-pion decay mode. This is because that in the single pion decay process, the leading 
non-vanishing contribution comes from the sub-leading diagrams shown in Fig. 1, where three vertices exist, whereas 
in the double-pion decay process, the leading non-vanishing contribution comes from the two vertices diagrams shown 
in Refs. [20, Bil. Clearly, the obtained very small decay width for the single pion decay mode is consistent with the 
current experimental status that no d* + N Nv process has been found in the present data set. 


Some approximations in calculation should be further discussed. In the sub-leading diagrams shown in Fig. 1, 
the Goldstone-Boson exchange between two nucleons must be introduced to convert two As to two nucleons. From 
PDG Bi, one finds that the largest decay mode for A(1232)3/2* is Na with a branching ratio of about 100%. 
Therefore, in the realistic calculation, considering the pion-exchange only would not miss the major feature and 
make the result meaningless. Due to the larger mass of nucleon excitations, we do not take them as the intermediate 
nucleon state. We also ignore the contribution from the diagrams where the intermediate A state exist, because the 
quark model calculation tells us that the AAz coupling faa« is much smaller than the N Nz coupling fNNr [23]. 
Again, we would specially emphasize that the contributions from the large CC component in the d* wave function 
could be ignored. The reason is the following. In the previous paragraph, we have mentioned the extracted channel 
wave functions for various channels have already absorbed the effect of the totally anti-symmetrization, and the 
channel wave functions for the AA and CC channels are orthogonal to each other. So in our decay calculation, 
the inter-cluster quark exchange should not be considered anymore. In the (SJ = 30) case, where S and I denote 
the spin and isospin of the system, respectively, converting CC to AA requires an exchange of a colored object, 
namely a gluon. The calculation shows that without quark exchange, the matrix elements of the one-gluon-exchange 
interaction (OGE) between AA and CC are zero, namely CC cannot be converted to AA. If one would convert 
CC(ST = 30) to AA(S 4 3) and CC to NN, because the spin of the system is S = 3 and the parity is positive, it 
needs at least D-wave between A and A with (S 4 3) and between N and N, respectively. Then, the required tensor 
force in OGE, which is a higher order term, would make these conversions suppressed greatly. In short, we can ignore 
the contribution from the CC component in d* in the d* decay calculation. 
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To summarize, we proceed a calculation for the single-pion decay mode of d* with the help of our wave function 
obtained in the chiral constituent quark model calculation. It shows in our calculations that the CC 
component has a large fraction of about 2/3 in the wave function of d*, and the rest part, the AA component, 
is responsible for its widths in the decays of d* — daa and d* — NWNzam as well as d* + NNr. As a result, 
the obtained partial widths for the double-pion decay modes of d* are in good agreement with the experimental 
measurement. Moreover, the obtained width for the single-pion decay model in this calculation is much smaller 
than those in the double-pion decay modes. If we assume that the total width of the d* is about 75 MeV, the 
predicted branching ratio of the single-pion decay mode d* — pnn? is about 3.0% in the one-channel case and 
0.9% in the coupled-channel case. It should be emphasized that the obtained single-pion decay width is much 
smaller than the widths of the double-pion decay modes. This result agrees with the present experiment status 
that such a single-pion decay mode has not been found in that data analysis. It is quite different from the 
width reported by the investigation with the AN7 assumption for the structure of d*, where the predicted width 
for the single pion decay mode would be large enough to be observed in the experimental measurement. It is ex- 
pected that an intensive data analysis of the d* + pnr? channel would judge different explanations for the nature of d*. 
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